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j M than the least-squares method when the results for the parameters were 
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; f S assessed using the goodness of fit measures. Also, we obtained the 
estimation, Weibull distribution. a ast od l o. 
asymptotic distribution of the MLE which was asymptotically efficient as 


the sample size increases. The inverse of the Fisher’s information matrix 
which is the asymptotic variance-covariance matrix was also obtained. 





I. INTRODUCTION (Razali et.al, 2009). Pobocikova et.al (2017) investigated 
The Weibull distribution has gained much weight in the performance of six estimators using the Monte Carlo 


the real world and is mostly used in reliability and lifetime simulation for wind speed, the result of the study indicates 


analysis or survival analysis. The 2-parameter Weibull that the maximum likelihood estimation method performs 


distribution has the shape parameter ( 8 ) and the scale better in estimating the parameters of the 2-parameter 


= Weibull distribution for wind speed when the bias and root 
parameter (7 or @ ) whose values affect the characteristic 
. a . . mean square error (RMSE) were compared. However, 
life of the distribution, the reliability function, and the according to Chu & Ke (2012), the least-squares method 


failure rate. Most distributions such as normal distribution, outwits the maximum likelihood approach when the 


gamma distribution, inverse gamma, and some other sample size is small. 


common distributions have two parameters which makes Th; i e3 Weibull 
the 2-parameter Weibull distribution of immense interest. B A A ee ae A a 
According to Sun(1997), the Weibull distribution does not 


belong to the 2-parameter exponential family neither are 


distribution by applying the Monte Carlo simulation to 
generate artificial data of failure times of a product 


a through which we obtained the parameter estimates of the 
the two parameters orthogonal. In estimating the uga Ww ý P 


parameters. OF ‘the, Weibull: bmo he method: of shape and scale parameters of the Weibull distribution by 


i h i likelih imati MLE 
moment achieved the best result when it was compared comparing the maximum likelihood estimation ( ) 


. l i h h f the li ion least- 
with other estimators such as the MLE and the LS using EO 200 EE OF NS neat Tepreeston Aca evans 


the mean square error (MSE) and the total deviation (TD) method. The two approaches of estimation were also 


fom The aparani ant: 7 parameter Weibull disnbiion applied to two different sets of data in reliability analysis. 
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The research also assesses the asymptotic efficiency of the parameters of the 2-parameter Weibull distribution. 
maximum likelihood estimator of the shape and scale 


II. MAXIMUM LIKELIHOOD ESTIMATOR 


The density function of the 2-parameter Weibull distribution is given as 





= p 
f (xX)= — exp (2) x>0,B>0,7>0 (1) 
7] 


The cumulative density function is given by 


F(x) =1-exp4- = (2) 
7 


p 
where £ is the shape parameter and 77 is the scale parameter. The reliability function is given by R(x) = ap- = | , the 
7] 


cumulative density function can be expressed in terms of the reliability function as F(x)=1- R (x) . 


Maximum likelihood estimation is a method used to estimate the parameter of a probability distribution. This is 


obtained by maximizing the likelihood function of the distribution. The maximum likelihood estimator of the Weibull 
distribution is given as follows: 





n n p-1 p 
TORTGAN f2) 
i=] i=] 


7] 


Mon B 
E) JI p-1 2) 
=|—; x: eXpy—| — (3) 
E i=l 1 


Taking the log of the likelihood in (3) we obtain the equation 


L(f(x)) = log P-nBlogiy +(B-1) Slows) | * | (4) 


i=] i=l 
The scale parameter can easily be obtained with an analytical method; the shape parameter / is very difficult to be obtained 


by the analytical method, however, it can be obtained by using the Newton-Raphson method or the least square method 
(Guureet al., 2012; Al-Fawzan, 2000). We obtain the scale parameter by differentiating (4) and solving the equation 


OL( f(x) _ 
on 
AL( f(x) _ BBS) ac o 
On 7 MK 
1 
yal | 
yn =| i=! 
n 
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Asymptotic efficiency 


The asymptotic efficiency of the 2-parameter Weibull distribution can be obtained as follows, We find the partial 
derivatives of the parameters by differentiating the log of the Weibull distribution. 


oy, -ÎLE SO -2 AED =) 
an? n nn n 


B 
olog f(x) 1 x 
ee ee | l 
a2 pP P i eef 
B B 
—Olog f(x) 1 {= A(z = 
= —*—— =——+—| >| +—]—]} log] — 
ev On Op n n\n n\n He 


» = Heel L,e 8) gl E)E] 
*! aB an n n\n n) n\n 


We obtain the determinant using the Hessian matrix of the Weibull distribution and ignoring all negative terms. 
H ( 0) = i O12 
CDi. 222 


where 0 = (2.7) 


p 2p 2 
1 +1 
|H (0) = cane + P 2 log 2 , considering log x < x+1 .We rewrite the determinant as 
2 2 
pn n n n 


7] 
moj LEAS ea] m 
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Considering = “10 f 


p 2p 
We obtain the upper bound of h Bi (x) from (7), for any x >0, =| < > + ) 
T No 


h 4(2B, (28,41 2h, 
) 2) A 2) (Za) (8) 
Pon 10 0 1o 
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Then E (h 5 
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B, 1, 


The kth moment is given as follows: 
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Therefore, based on Bahadur (1964) we have 


Vn(8,-8,)—>n, (01, OT) 


n n 


Hence we obtain the Fisher’s information matrix J, (0 ) and its inverse, the information matrix is given as 
(0) --2( E 
r 22 

For us to obtain the Fisher information in a simplified form, we make the same assumption made by Sun (1997) that for 721, 

2 , p 
we define y; = f [log( y)| exp(— y)dy , and we let Y = = be an exponential random variable with mean 1 whose density 

7] 

0 

function is expressed as f(y)=exp(—y), y20 


BE Bi loe[ž] = EZOO 


0 


([1+os(»)]exp(-y)as =1+y, 


rE BICAJ - [ris Fo) 
[feo +2log y Jexp(-y)dy = y3 +27) 


The Fisher’s information matrix can be expressed as 


Bo bn 
2 
(a)=| ' ! 
n p’ 
y2- 
The determinant of the information matrix pi (9) = a, where —y, is Euler’s constant and y, — Ve is the variance of 


7] 


log (g) where @ is an exponential random variable with a mean of 1. 





2 
[7(9) | = 4 ‘ 
yy-n | 1+7% po 
7] n’ 


II. LEAST-SQUARES ESTIMATION 
The least-squares estimation of the 2-parameter Weibull distribution is also important and can be used to estimate the 
parameters of the Weibull distribution inasmuch as the maximum likelihood is used. In estimating the parameter of the 
Weibull distribution, we will consider the cumulative density function and obtain a regression equation from which we can 
deduce the equations for estimating the parameters of the 2-parameter Weibull distribution. The cumulative density function 
of the Weibull distribution in (2) becomes: 
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ap- (2) J=1-Fe 


Taking the log of the reliability function, we obtain 


~ =—In[I-F(x)] 
n 


nA 


Sie 


From (9), let Y, =In(x;), X =In{-In [1— F(x)]} 


—In[1- F(x)]} 


In{—In[1— F(x)]! + In(77) (9) 


1 
a= 3 and b=In(7), we write (9) as a linear function of the form 


Y; =aX +b. By substituting Y, and X into the least-squares formula, we will obtain the estimates for a and b through which 


we estimate the parameters p and y. 


(10) 


(11) 


1 
hence we have p =— and 7 = exp(b) , where F (x) can be estimated by the mean rank approach (Pobocikova,2012), 
a 


IV. SIMULATION AND ANALYSIS 

In this study, we first simulated data using the Monte 
Carlo simulation. Random datasets were generated from 
the uniform distribution for the study. The datasets 
generated from the uniform distribution were sampled 
uniformly from the range of O and 1 to represent the 
failure times(hours) for a particular product. The randomly 
selected samples from the uniform distribution were 


simulated with the model X = n(—In (<) A to represent 


the failure times for the product as a random sample from 
the Weibull distribution, where z has the uniform 
distribution (0,1). We simulated 36,000 failure times and 
chose n = 50,100 , and 1000, we chose different scale and 


shape parameters to simulate the data for the Weibull 
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distribution, and the sample sizes n were chosen after 
each simulation, this was replicated 1000 times in each 
simulation. The shape parameter values G were chosen so 


that 6 =0.5 will represent the scenario when (<1 and 
the other values (2,3) of the shape parameter when £ >1, 
and £=1 when we have constant failure times. The two 


parameters for the Weibull distribution were estimated 
through the use of the maximum likelihood estimation and 
the least-squares method. The performance of these two 
estimation techniques was compared using the goodness of 
fit measures. The maximum likelihood approach estimates 
the parameters of the 2-parameter Weibull distribution 
better than the least-squares method. The results for the 
parameter estimation using the maximum likelihood and 
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the least-squares method are tabulated in Table 1. By 
comparing the goodness of fit for the MLE and the LS, the 


the Weibull distribution for this study are shown in Fig.1. 
The probability plot and the cumulative density function 


MLE performs better in estimating the model parameters 


than the least-squares method. 


The probability density function and the cumulative 


density of the different parameter values used to simulate 


(CDF) are shown in Fig. 2. 


Table 1. Parameter estimation of the Weibull distribution 


Maximum Likelihood Estimation 


Parameter 
Values Estimates Confidence Interval Goodness of fit 
1) Log- 
ae h AA h BI AD 
je r 2 1 (œ (SE) B (SE) Cy (a )) CK £) likelihood  AICe £ 
50 2 0.5 1.8310(0.5454) 0.4949(0.0544)  (1.0213,3.2829) (0.3990,0.6138) -87.4434 179.142 182.711 0.4928 
1 1 1.2068(0.1495) 1.2060(0.1315)  (0.9467,1.5384) (0.9739,1.4934) -54.8121 113.880 117.448 0.4579 
3 2 3.2545(0.2276)  2.1122(0.2389)  (2.8377,3.7325) (1.6923,2.6363) -88.0339 180.323 183.892 1.0994 
2 3 2.0121(0.0952) 3.1467(0.3563) (1.8340,2.2075) (2.5205,3.9285) -47.8072 99.8696 103.438 0.4987 
100 2 0.5 1.9099(0.4687)  0.4288(0.0344)  (1.1807,3.0896) (0.3664,0.5018) -170.569 345.261 350.348 0.5962 
1 1 1.2241(0.1171) 1.1052(0.0853)  (1.0148,1.4766) (0.9501,1.2856) -115.496 235.115 240.202 0.7073 
3 2 2.8691(0.1373)  2.1891(0.1756) (2.6122,3.1512) (1.8706,2.5617) -160.233 324.589 329.676 0.8453 
2 3 1.9377(0.0720)  2.8296(0.2240)  (1.8017,2.0841) (2.4229,3.3045) -100.351 204.826 209.912 0.3510 
1000 2 0.5 1.9012(0.1252)  0.5065(0.0125)  (1.6710,2.1631) (0.4826 ,0.5315) -1767.24 3538.5 3548.3 0.7730 
1 1 1.0413(0.0345) 1.0042(0.0247)  (0.9757,1.1112) (0.9569,1.0538) -1038.67 2081.36 2091.16 0.2608 
3 2 3.0284(0.0495)  2.0362(0.0505)  (2.9329,3.1270) (1.9395,2.1376) -1692.99 3389.98 3399.79 0.2070 
2 3 = 2.0028(0.0212) 3.1398(0.0776) (1.9616,2.0449) (2.9913,3.2956) -944.250 1892.51 1902.32 0.3783 
Least-squares Estimation 
7] Log- 
ara À AnA D BI AD 
we (Q) £ 1 (@ ) (SE) P (SE) Cl(77 (@ )) Cl( B ) likelihood  AJCc . 
50 2 0.5 1.8310(0.5454) 0.4949(0.0544)  (1.0213,3.2829) (0.3990,0.6138) -87.4434 179.142 182.711 0.4928 
1 1 1.2010(0.1463) 1.2223(0.1329)  (0.9459,1.5248) (0.9877,1.5125) -54.8233 113.902 117.471 0.4810 
3 2 3.3544(0.3018) 1.7234(0.2069) (2.8121,4.0014) (1.3620,2.1806) -89.7926 183.84 187.409 1.5810 
2 3 =. 2.0162(0.1016) 2.9541(0.3382)  (1.8266,2.2254) (2.3603,3.6972) -47.9805 100.216 103.785 0.4353 
100 2 0.5 1.9055(0.4938)  0.4123(0.0333)  (1.1467,3.1666) (0.3519,0.4831) -170.697 345.518 350.604 0.4999 
1 1 1.2004(0.1067) 1.1784(0.0897)  (1.0085,1.4289) (1.0150,1.3680) -116.000 236.124 241.21 1.0681 
3 2 2.9312(0.1672) 1.8901(0.1571)  (2.6211,3.2779) (1.6059,2.2245) -162.173 328.469 333.556 1.1726 
2 3 1.9515(0.0784)  2.6302(0.2117)  (1.8037,2.1114) (2.2463,3.0798) -100.866 205.855 210.941 0.3884 
1000 2 0.5 1.8555(0.1157) 0.5229(0.0127)  (1.6420,2.0967) (0.4985 ,0.5485) -1768.47 3540.94 3550.75 1.4114 
1 1 1.0383(0.0341) 1.0122(0.0249)  (0.9736,1.1076) (0.9646,1.0621) -1038.74 2081.5 2091.31 0.3034 
3 2  3.0364(0.0507) 1.9972(0.0498)  (2.9387,3.1374) (1.9019,2.0972) -1693.37 3390.76 3400.56 0.2147 
2 3 =. 2.0024(0.0213)  3.1323(0.0775) (1.9611,2.0445) (2.9841,3.2879) -944.255 1892.52 1902.32 0.3718 
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Weibull Distribution 





IA: Probability density plots for the different parameter values 


Weibull Distribution 


— ßf=0.5, n=1 
--- B=1.0, n7=1 





1B: Cumulative density plots for the different parameter values 


Fig.l: Plots of the probability density and cumulative density of the Weibull distribution with different values of the shape 
(p ) and scale parameter (n (œ )). 


50 samples 100 samples 1000 samples 
99.9% A 


99% 
95% 
90% 
80% 
70% 
60% 
50% 
40% 


30% 


20% 


Fraction failing 
Fraction failing 
Fraction failing 


10% 


10% 


5% 


3% 


2% 


--- True CDF --- True CDF 


1% — Fitted Weibull_2P (a=2.849, B=1.757) ~ 0.5% — Fitted Weibull _2P (a=2.81, B=1.892) 


=-=- True CDF 
— Fitted Weibull _2P (a=2.934, B=2.036) 





Time 


Fig.2: Plot of the cumulative density function and the probability plot of the data values of three different samples randomly 
selected from the Weibull distribution and their estimated values of the shape ( p ) and scale parameter (n(@ )). 
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4.1 Applications in Reliability Analysis 
Health 


Considering the 
distribution in reliability analysis, we use the data on 128 
randomly selected remission times (months) of bladder 


application of Weibull 


cancer patients. The data have been used by many 
researchers for different research activities, Zea et al. 
(2012) used the data comparing sub-model including beta- 
Pareto distribution and others, and the beta exponential 
Pareto distribution. Almheidat et al. (2015) also 
investigated the Cauchy-Weibull distribution with this 
dataset. We apply the Weibull distribution to fit the same 
data, however, we removed all outliers that are 2-standard 
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deviation away from the mean from the dataset before 
fitting the Weibull distribution to the data. In order to 
ascertain that the correct distribution is fitted, we fitted 
many other distributions such as the Lognormal, Gamma, 
Exponential, Loglogistic, etc. to the data after removing all 
outliers from the data and then select the best distribution 
that fits the data. This reduces the sample size from 128 to 
122 remission times of bladder cancer patients. Fig.3 
indicates the fit of distributions to the data. The data for 
the 128 remission times (months) of bladder cancer 
patients are shown below in Table 2. 


Table 2: Remission times (months) of bladder cancer patients 


0.08 1.46 2.69 3.57 
11.98 15.96 25.82 0.2 
5.71 7.39 9.02 12.02 
2.15 3.7 4.98 5.85 
32.15 0.5 2.02 2.83 
9.47 12.07 17.14 34.26 
5.09 6.54 7.63 9.74 
2.09 3.02 4.18 5.17 
18.1 43.01 0.9 PAPIS 
7.87 10.34 13.29 19.13 
4.26 J2 6.94 7.93 
1.19 2.46 3.36 4.33 
14.24 21,13 1.26 2.54 
8.37 11.25 14.76 22.69 
5.41 7.26 8.53 11.64 
3.32 4.5 5.49 7.28 
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4.51 5.62 da2 8.66 
1.76 2.69 3.64 4.87 
16.62 26.31 0.4 2.02 
7.59 9.22 12.03 17.12 
3.82 5.06 6.25 7.62 
0.51 2.07 2.87 3.88 
12.63 17.36 36.66 0.81 
6.76 7.66 10.06 13.11 
3.25 4.23 5.32 6.93 
46.12 1.05 2.26 3.31 
10.66 13.8 20.28 79.05 
5.34 6.97 8.26 10.75 
3.36 4.34 5.41 7.09 
1.35 2.62 3.48 4.4 
14.77 23.63 1.4 2.64 
8.65 11.79 14.83 25.74 
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Histogram plot of each fitted distribution 


Probability Density Function 





Data 


Weibull (a, 8) 
Gamma (a, B) 
Weibull (a, B, Y) 
Gamma (a, B, Y) 

— Exponential (A) 
Exponential (A, Y) 
Loglogistic (a, 8) 
Loglogistic (a, B, y) 
Lognormal (u, a) 


Lognormal (uz, o, Y) 


Normal (u, g) 
Gumbel (u, 7) 


Cumulative probability density 


Cumulative Distribution Function 





Data 


Fig.3: Histogram with each fitted distribution with the plots of its PDF and CDF. 


Distribution Alpha 

Weibull 2P 8.204 
Gamma_2P 5.130 
Weibull 3P 8.128 
Gamma_3P 5.130 


Exponential_1P 
Exponential_2P 
Loglogistic_2P 
Loglogistic_3P 
Lognormal_2P 
Lognormal_3P 
Normal_2P 
Gumbel_2P 


Table 3: Comparison of the various distributions for the remission times data 


Beta 


1.276 
1.482 
1.258 
1.482 


1.864 
1.864 


Gamma 


0.0413 


Mu Sigma 
1.655 0.998 
1.655 0.998 
7.603 6.021 

10.922 7.362 


Lambda 


0.13153 
0.13293 


Log- 
likelihood 


-364.19 
-364.34 
-364.11 
-364.34 
-369.77 
-368.19 
-370.61 
-370.61 
-374.72 
-374.72 
-392.13 
-420.56 


AICc 


732.481 
732.787 
734.421 
734.889 

741.58 
740.483 
745.312 
747.414 
753.539 
755.642 

788.35 
845.226 


Weibull (a, 8) 
Gamma (a, 8) 
Weibull (a, 8, Y) 
Gamma (a, B, Y) 
Exponential (A) 
Exponential (A, y) 
Loglogistic (a, B) 
Loglogistic (a, B. Y) 
Lognormal (u, 0) 


Lognormal (4, g, Y) 


Normal (u, 0) 


Gumbel (u, 7) 





BIC 


737.988 
738.294 

742.63 
743.098 
744.351 

745.99 
750.819 
755.623 
759.046 

763.85 
793.857 
850.734 


AD 


0.322 
0.262 
0.315 
0.262 
3.479 

2.16 
0.566 
0.566 
1.365 
1.365 
4.218 

8.34 


Table 3 compares the various distributions using the AICc, and BIC and it was observed that the 2-parameter 


Weibull distribution is the best fit for the data. Fig. 4 shows the probability plot of each of the distributions compared, from 


Fig.4 there is an indication that the 2-parameter Weibull distribution fits better than the other distributions. 
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Probability plots of each fitted distribution 








Weibull_2P Gamma_2P Weibull_3P Gamma_3P 
Exponential_1P Exponential _2P Loglogistic_2P Loglogistic_3P 
Lognormal_2P Lognormal_3P Normal_2P Gumbel_2P 





Fig.4: Probability plot of the remission times (months) of the bladder cancer patients for the various distributions. 


The 2-parameter Weibull distribution is the Weibull distribution, the results of the two methods are 


considered as the best distribution for the data, we, 
therefore, study the reliability of the remission times of the 
bladder cancer patients using the distribution chosen. The 


compared using the AICc and the BIC. The parameter 
estimates are shown in Table 4. The MLE approach 
performs better in estimating the parameters of the 


maximum likelihood estimation method and the least- remission times of the bladder cancer patients compared 


squares method are considered in fitting the parameters of with the least-squares method. 


Table 4: Parameter estimate of the Weibull distribution for the remission times 


Point Standard Log- 
Parameter Estimate Error Lower CI Upper CI likelihood AICc BIC AD 
MLE Alpha 8.20364 0.613436 7.08528 9.49852 -364.19 732.481 737.988 0.322 
Beta 1.27557 0.0898488 1.11108 1.4644 
LS Alpha 8.24435 0.63076 7.0963 9.57813 -364.24 732.576 738.084 0.317 
Beta 1.25165 0.08858 1.08953 1.43789 
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Probability plot 
Weibull CDF 


— Fitted Weibull 2P (a=8.204, B=1.276) 


Fraction failing 


10% 


1% 


01 1 10 
Time 


Fraction surviving 
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Weibull Distribution 
Survival Function 
a=7,B=1.0 


— fitted distribution 
=== original distribution 


10 


0.4 


co 
mi 


04 


0z 





0.0 F 
0 10 20 30 40 


x Values 


Fig.5: Probability plot and the plot of the reliability function for the remission times. 


The plot in Fig.5 gives the probability plot for 
the Weibull distribution and the reliability function for the 
remission times (month) for the bladder cancer patients. 
The estimated shape parameter (1.276) using the MLE 
indicates that the remission rate of the bladder cancer 
patient increase as the remission time increases. The B10 
life of the data was 1.4 months which indicates that we 
expect 10% of the bladder cancer patients to be remitted 
after 1.4 months. 


Industrial 


We again consider breaking stress of carbon 
fibers of 50mm length (GPa), the data was first used by 


Nicholas and Padgett (2006), Lemonte (2011) using £ - 


Birnbaum-—Saunders distribution, Al-Agtash et.al (2014) 
use the data to study the application of Gumbel-Weibull 
distribution. For our study, we remove all outliers that are 
2- standard deviation away from the mean, this reduces the 
dataset from 66 to 59. The dataset was fitted with other 
distributions and the best distribution was the 2-parameter 
Weibull distribution. Fig.6 illustrates the fit of the 
distributions to the breaking stress of the carbon fibers data. 
The carbon fibers data is shown in Table 5. 


0.39 
1.87 
2.53 
2.85 
3.15 
3.56 
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0.85 
1.89 
239 
2.87 
3.19 

3.6 


1.08 
2.03 
299 
2.88 
Be 
3.65 


Table 5: Breaking stress of Carbon fibers 


1.25 
2.03 
2.56 
2.93 
3.22 
3.68 


1.47 
2.05 
2.59 
2.95 
3.27 

N, 


1.57 1.61 1.61 1.69 1.8 1.84 
212 2.35 2.41 2.43 2.48 2.5 
2.67 2.73 244 2419 2.81 2.82 
2.96 297 3.09 3.11 3.11 3.15 
3.28 3.31 3.31 3.33 3.39 3.39 
210 4.2 4.38 4.42 4.7 4.9 
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Probability Density Function 


— Weibull (a, 8) 


Normal (u, 7) 
Gumbel (u, g) 


Weibull (a, 8, Y) 
— Gamma (a, B) 
Lognormal (u, 0) 
Gamma (a, B, Y) 
Loglogistic (a, B) 
Lognormal (u, g, Y) 


Loglogistic (a, B, Y) 
Exponential (A, y) 


Exponential (A) 


Histogram plot of each fitted distribution 


Cumulative Distribution Function 
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Normal (u, 0) 
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Weibull (a, 8, Y) 


© 
n 


Gamma (a, B) 
Lognormal (u, 0) 
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Loglogistic (a, B, y) 
Exponential (A, Y) 
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0.2 





Fig.6: Histogram with each fitted distribution with the plots of its PDF and CDF. 


Exponential (A) 





Table 6: Comparison of the various distributions for the breaking stress of carbon fibers data 


auto Alpha Beta Gamma Mu Sigma 
Weibull_2P 2.990 4.836 

Normal_2P 2.736 0.661 
Gumbel_2P 3.056 0.591 
Weibull_3P 2.990 4.836 0.000 

Gamma_2P 0.180 15.203 

Lognormal_2P 0.973 0.268 
Gamma_3P 0.180 15.203 0.000 

Loglogistic_2P pS pe, 6.619 

Lognormal_3P 0.000 0.973 0.268 
Loglogistic_3P 2422 6.619 0.000 

Exponential_2P 1.250 


Exponential_1P 


Table 6 compares the distributions and the best 
distribution that best fit the breaking stress of the carbon 
fibers data after removing all observations that were 2- 
standard deviation away from the mean was the 2- 
parameter Weibull distribution. The comparison was done 
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Log- 
likelihood 


-58.2972 
-59.3269 
-59.9252 
-58.2972 
-61.4958 
-63.3661 
-61.4958 
-63.6408 
-63.3661 
-63.6408 
-82.3637 
-118.378 


AICc 


120.809 
122.868 
124.065 
123.031 
127.206 
130.946 
129.428 
131.496 
133.169 
133.718 
168.942 
238.827 


BIC 

124.749 
126.809 
128.005 
128.827 
131.147 
134.887 
135.224 
135.437 
138.965 
139.514 
172.883 
240.834 


AD 
0.578 
0.771 
0.519 
0.578 
1.265 
1.595 
1.265 
1.192 
1.595 
1.192 
7.952 

15.551 


by using the AICc and BIC. From Fig.7, it can be observed 
that the probability plot of the 2-parameter Weibull 


distribution achieves a better result than the other 


distributions. 
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Probability plots of each fitted distribution 


Weibull_ 2P Normal_2P 


Gumbel_2P Weibull_ 3P 





Gamma_2P Lognormal_2P 





Lognormal_3P Loglogistic_3P 





Gamma_3P Loglogistic_2P 





Exponential_1P 


Exponential_2P 


Fig.7: Probability plot of breaking stress of carbon fibers data for the various distributions. 


We study the breaking stress of carbon fibers 
using the Weibull distribution to assess its reliability. We 
use the maximum likelihood estimation and the least- 
squares method to estimate the scale and shape parameter 


of the carbon fibers data, we compare these two methods 
using the AICc and the BIC. Table 7 presents the 
parameter estimate of the breaking stress of the carbon 
fibers using the 2-parameter Weibull distribution. 


Table 7: Parameter estimate of the Weibull distribution for the breaking stress of carbon fibers 


Point Standard 


Parameter re Tan Lower CI 
MLE Alpha 2.99031 0.08466 2.8289 
Beta 4.83619 0.50136 3.94694 
LS Alpha 2.99167 0.08904 2.82214 
Beta 4.83619 0.50136 3.94694 
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Log- 

Upper Cl ican , AlCe BIC AD 
3.16093 -58.297 120.809 124.749 0.578 
5.92579 
3.17138 -58.436 121.086 125.027 0.547 
5.92579 
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Fig.8: Probability plot and the plot of the reliability function for the breaking stress of carbon fibers. 


The plot in Fig.8 gives the probability plot for the Weibull 
distribution and the reliability function for the breaking 
stress of the carbon fibers data. The maximum likelihood 
estimator gives the best estimate with an estimated shape 
parameter (4.836) which indicates that the rate of breaking 
stress increase as the stress time increases. The B10 life of 
the data was 2.77206 GPa which indicates that we expect 
10% of the carbon fibers to break after obtaining breaking 
stress of 2.77206 GPa. 


V. CONCLUSION 

The paper examines two approaches for estimating the 
parameters of the 2- parameter Weibull distribution. In this 
research, we compared two estimation methods using 
artificial data generated with Monte Carlo simulation and 
two real datasets. The research findings indicate that the 
maximum likelihood estimation method performs better 
than the least-squares approach. 


www.ijaers.com 


In general, the Weibull distribution is very useful in 
many areas such as survival analysis, reliability analysis, 
etc. The estimation of the 2-parameter or 3-parameter 
Weibull distribution with or without censored data has 
gained much attention in these areas of study. In this paper, 
we compared the maximum likelihood estimation and the 
least-squares method, and by assessing the properties of 
the distribution and the estimates of the 2-parameter 
Weibull distribution using the goodness of fit measures. 
We observed that the maximum likelihood estimation 
performs better than the least-squares method and the ML 
estimator is asymptotically efficient as it converges in 
distribution as the sample size increases. 


In order to investigate the application of the Weibull 
distribution in reliability analysis. We estimated the 
parameter of the Weibull distribution considering two 
different datasets using the MLE and LS and examined the 
failure rate of the datasets. The findings also indicate that 
the MLE performs better than the least-squares method. 


Page | 223 


Edward Appau Nketiah et al. International Journal of Advanced Engineering Research and Science, 8(9)-2021 


REFERENCES 


[1] Al-Fawzan, M. A. (2000). Methods for estimating the 
parameters of the Weibull distribution. King Abdulaziz City 
for Science and Technology, Saudi Arabia. 

[2] Al-Agtash, R., Lee, C., & Famoye, F. (2014). Gumbel- 
Weibull distribution: Properties and applications. Journal 
of Modern applied statistical methods, 13(2), 11. 

[3] Almheidat, M., Famoye, F., & Lee, C. (2015). Some 
generalized families of Weibull distribution: Properties and 
applications. International Journal of Statistics and 
Probability, 4(3), 18. 

[4] Bahadur, R. R. (1964). On Fisher's bound for asymptotic 
variances. The Annals of Mathematical Statistics, 35(4), 
1545-1552. 

[5] Cordeiro, G. M., & Lemonte, A. J. (2011). The B- 
Birnbaum-—Saunders distribution: An improved distribution 
for fatigue life modeling. Computational Statistics & Data 
Analysis, 55(3), 1445-1461. 

[6] Chu, Y. K., & Ke, J. C. (2012). Computation approaches 
for parameter estimation of Weibull 
distribution. Mathematical and computational 
applications, 17(1), 39-47. 

[7] Guure, C. B., Ibrahim, N. A., & Ahmed, A. O. M. (2012). 
Bayesian estimation of 2-parameter weibull distribution 
using extension of Jeffreys' prior information with three 
loss functions. Mathematical Problems in 
Engineering, 2012. 

[8] Nicholas, M. D. & Padgett, W. J. (2006). A bootstrap 
control for Weibull percentiles. Quality and Reliability 
Engineering International, 22, 141-151. 

[9] Pobocikova, I., & Sedliackova, Z. (2012). The least square 
and the weighted least square methods for estimating the 
Weibull distribution parameters—a comparative 
study. Communications-Scientific letters of the University 
of Zilina, 14(4), 88-93. 

[10] Pobocikova, I., Sedliackova, Z., Michalkova, M., & George, 
F. (2017). Monte Carlo Comparison of the Methods for 
Estimating the Weibull Distribution Parameters—Wind 
Speed Application. Communications-Scientific letters of 
the University of Zilina, 19(2A), 79-86. 

[11] Razali, A. M., Salih, A. A., & Mahdi, A. A. (2009). 


Estimation accuracy of Weibull distribution 
parameters. Journal of Applied Sciences Research, 5(7), 
790-795. 


[12] Sun, D. (1997). A note on noninformative priors for 
Weibull distributions. Journal of Statistical Planning and 
Inference, 61 (2), 319-338. 

[13] Zea, L. M., Silva, R. B., Bourguignon, M., Santos, A. M., 
& Cordeiro, G. M. (2012). The beta exponentiated Pareto 
distribution with application to bladder cancer 
susceptibility. International Journal of Statistics and 
Probability, 1(2), 8. 


www.ijaers.com Page | 224 


